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ABSTRACT 

Weak gravitational lensing observations probe the spectrum and evolution of density fluctuations 
and the cosmological parameters which govern them. The non-linear evolution of large scale structure 
produces a non-Gaussian signal which is potentially observable in galaxy shear data. We study the three- 
point statistics of the convergence, specifically the bispectrum, using the dark matter halo approach which 
describes the density field in terms of correlations between and within dark matter halos. Our approach 
allows us to study the effect of the mass distribution in observed fields, in particular the bias induced 
by the lack of rare massive halos (clusters) in observed fields. We show the convergence skewness is 
primarily due to rare and massive dark matter halos with skewness converging to its mean value only if 
halos of mass M > 10 15 Mq are present. This calculational method can in principle be used to correct for 
such a bias as well as to search for more robust statistics related to the two and three point correlations. 

Subject headings: cosmology: theory — large scale structure of universe — gravitational lensing 
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et al 1997) or numerical ly through simulations (e.g., Jain 



We ak gravitational lonoing of faint galaxioo proboo the 

distribution of matter along the line of sight. Lensing 
by large-scale structure (LSS) induces corr elation in the 
galaxy e llipticities at the percent level (e.g., Blandford ct 
g 1991; |Miralda-Escude] 1991; [Kaiscr| 1992). Though chal- 
lenging to measure, these correlations provide important 
cosmological information that is complementary to that 
supplied by the cosmi c microwave ba ckgro und and poten- 
tially as precise (e.g., Jain fc Seljak 1997; Bernardcau et 
|aj 199 7; |Kaiser| 1998 ; [Schneider et aj 199 8; |Hu fc Tegm~ 
1999; ICooravl 1999 ; [Van Waerbeke et al| 1999; see [Bartel- 
|mann fc Schneider 2000 for a recent review). Indeed sev- 
eral recent studies have provided the first cle ar evidence 
for weak le nsing in so-called b lank fi elds (e.g., IVan Waer 



beke ct all 2 000; [Bacon et a| 2000; |Wittman ct al| 2000; 



Kaiser et al 2000), though more work is clea rly needed to 
understand even the statistical errors (e.g. Cooray et aj 



2000b). 

Given that weak gravitational lensing results from the 
projected mass distribution, the statistical properties of 
weak lensing convergence reflect those of the dark mat- 
ter. Non-linearities in the mass distribution induce non- 
Gaussianity in the convergence distribution. With the 
growing observational and theoretical interest in weak 
gravitational lensing, statistics such as the skewness have 
been suggested as probes of cosmological parameters and 
the non-linear e voluti on of larg e scal e st ructu re (e.g., 
Bcrnardeau et all 1997; [Jain et all 2000; |lu] 1999; |Munshj 



fc Jain| 1999j |Van Waerbeke et al| 1999) 



Here, we extend previous studies by considering the full 
convergence bispectrum, the Fourier space analog of three- 
point function. The bispectrum contains all the informa- 
tion present at the three point level, whereas conventional 
statistics, such as skewness, do not. The calculation of the 
convergence bispectrum requires detailed knowledge of the 
dark matter density bispectrum, which can b e obtain an 



alytically through perturbation theory (e.g., Bernardeau 

Alfred P. Sloan Fellow 



et al 



2000; 



White fc Hu| 1999). Perturbation theory, how- 



is not applicable at all scales of interest, while nu- 
merical simulations are limited by computational expense 
to a handful of realizations of cosmological models with 
modest dynamical range. Here, we use a new approach 
to obtain the density field bispectrum analytically by de- 
scribing the underlying three point correlations as due to 
contributions from (and correlations between) individual 
dark matter halos. 

Techniques for studying the dark matter density field 
through halo contributio ns hav e recently been dev eloped 
( |Seljak| 2000; |Ma fc Fry| 2000b; [Bcoccimarro ct alj 2000) 
and applied to two-point lensing statistics ( Gooray ct aj 
2000b). The cr itical ingredients a re: the Press-Schcchtcr 
formalism (PS; Press fc Schcchterj 197 4) for the mass func- 



tion; the NFW profile of Navarro et~a] (1996), and the halo 
bias model of Mo et al. (1997). The dark matter halo 



approach provides a physically motivated method to cal- 
culate the bispectrum. By calibrating the halo parameters 
with N-body simulations, it can be made accurate across 
the scales of interest. Since lensing probes scales rang- 
ing from linear to deeply non-linear, this is an important 
advantage over perturbation-theory calculations. 

Throughout this paper, we will take ACDM as our fidu- 
cial cosmology with parameters il c = 0.30 for the CDM 
density, flf, = 0.05 for the baryon density, fl\ = 0.65 for 
the cosmological constant, h = 0.65 for the dimensionless 
Hubble constant and a scale invariant spectrum of primor- 
dial fluctuatio ns, normalized to galaxy cluster abundances 
(o"8 = .9 see Viana fc Liddlc 1999) and consistent with 
COBE flBunn fc White| 1997). For the linear power spec- 



trum, w e take the fitting formula for the transfer function 
given in [Eiscustcin fc Hu| (1999). 

In we review the dark matter halo approach to mod- 
eling the density field. In §|| we apply the formalism to the 
convergence power spectrum, skewness, and bispectrum. 
We summarize our results in §||. 

2. DENSITY POWER SPECTRUM AND BISPECTRUM 
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2.1. General Definitions 

Underlying the halo approach is the assertion that dark 
matter halos of virial mass M are locally biased tracers of 
density perturbations in the linear regime. In this case, 
functional relationship between the overdensity of halos 
and mass can be expanded in a Taylor series 

5 h (x, M; z) = bi(M; z)S(x; z) + ^b 2 {M; z)6 2 (x; z) + . . . 

. " (1) 

Mo et al. (1997) give the following analytic predictions for 
the bias parameters which agree well with simulations: 



h.{M\z) = 1 + 



v 2 {M;z) - 1 



and 



b 2 (M;z) = -[h(M;z)-l] 



v 2 {M; z)-3 
<j 2 (M;z) 



(2) 



(3) 



Here v(M, z) — S c /a(M, z), where cr(M, z) is the rms fluc- 
tuation within a top-hat filter at the virial radius corre- 
sponding to mass M, a nd S c is the threshold overdensity of 
spherical collapse (see Henry 2000 for useful fitting func- 
tions). 

Roughly speaking, the perturbative aspect of the clus- 
tering of the dark matter is described by the correlations 
between halos, whereas the nonlinear aspect is described 
by the correlations within halos, i.e. the halo profiles. We 
will consider the Fourier analogies of the 2 and 3 point 
correlations of the density field defined in the usual way 



(4) 



(a(ki)a(k 2 )a(k3)) = (2nfs(k 1 + k 2 + k 3 )B\k u k 2 , k 3 ) . 

(5) 

Here and throughout, we occasionally suppress the redshift 
dependence where no confusion will arise. 

As we shall see, these spectra are related to the linear 
density power spectrum P{k) through the bias parame- 
ters and the normalized 3d Fourier transform of the halo 
density profile p(r, M) 



y(k,M) = j I 



dr 4nr 2 p(r, M) 



sin(fcr) 
kr 



(6) 



Note that y(k, M) can be written as a combination of 
sine and cosine integrals for computational purposes and 
y(k, M) -> as k -> 0. 

It is convenient then to define a general integral over the 
halo mass function dn/dM, 

I?{k l ,...,k ll ;z) = J dM^y ^L(M,z)bp(M) 

xyikuMj.-.yik^M), (7) 



where 6n = 1. We use the Press-Schechter (PS; Press & 



Schec iter 1974) mass function to describe dn/dM. Wc 
take the minimum mass to be 10 3 Mq while the maxi- 
mum mass is varied to study the effect of massive halos 
on lensing convergence statistics. In general, masses above 
10 16 Mq do not contribute to low order statistics due to 
the exponential decrease in the number density of such 
massive halos. 



2.2. Power Spectrum and Bispectrum 



Following Seljak (2000), we can decompose the density 
power spectrum, as a function of redshift, into contribu- 
tions from single halos (shot noise or "Poisson" contribu- 
tions), 

P pp (fc)=/ 2 °(M), (8) 
and correlations between two halos, 



p hh (fc) = [i\{k)yp{k) 



such that 



P t = P 



pp 



P 



hh 



(9) 
(10) 

As k 0, P hh -► P(k). 

Similarly, we decompose the bisp ectrum into terms in- 
volvi ng one, two and three halos (see Bcherrer & Bertschinger 
1991; [Via &; Fry) 2000b): 



B t =B 



ppp 



B 



Phh 



B 



hhh 



where 

B ppp (fci,*2,fc3) ^(fci^fcg). 
for single halo contributions, 



(11) 
(12) 



B phh (fc 1 ,fc 2 ,A: 3 ) =/ 2 1 (fc 1 ,fc 2 )/ 1 (fc 3 )F(fc 3 )+Perm. (13) 

for double halo contributions, and 

B hhh (ki, k 2 , k 3 ) = [2 J(h, k 2 , k 3 )ll(k 3 ) + I 2 (k 3 )] (14) 
x/ 1 1 (fc 1 )/ 1 (fc 2 )P(fc 1 )P(fc 2 ) + Perm. 

for triple halo contributions. Here the 2 permutations 
are k 3 <- > fci , k 2 . Second order p erturb ation theory tells 
us that ([Fry] 1984; |Bouchct ct~aj 1992; |Kamionkowski & 
Buchalter| 1999) 



j(k u k 2 ,k 3 ) = i--n m 2 / 63 + 



k 2 



k\ 



k 2 



k\ 



^2 



k 2 



k\ 



k 2 



2k x k 2 
'q-2/63 
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(15) 



As k — > 0, B hhh — > B PT , where B PT is the bispectrum 
predicted by second-order perturbation theory 



(16) 



B PT (fc!, k 2 ,k 3 ) = 2J{k 1 ,k 2 ,k 3 )P{k 1 )P{k 2 ) + Perm. 

with permutations following k 3 <-> k\, k 2 . 

2.3. Halo Profiles 

To apply the halo prescription for the power spectrum 
and bispectrum, we need to k now the profi le of the ha- 
los. We take the NFW profile ( Navarro et al 1996) with a 
density distribution 



p(r,M) 



{r/r s ){l + r/r s y 



(17) 



The density profile can be integrated and related to the 
total dark matter mass of the halo within r„ 



M = Airp s rl 



log(l 



(18) 



3 




k (h Mpc 



k(h Mpc ) 



Fig. 1. — Dark matter density spectra under the halo prescription, (a) Dark matter power spectrum at the present, (b) Mass cut off 
effects on the power spectrum, (c) Equilateral bispectrum at the present, (d) Mass cut off effects on the equilateral bispcctrum. The power 
spectrum is compared with the PD fitting function and the linear P(k) in (a), fn (b) and (d), from bottom to top, the maximum mass used 



in the calculation is 10 11 , 10 12 , 10 13 , 10 14 , 10 15 and fO 16 M 



0- 



where the concentration, c, is defined as r v /r s . Choos- 
ing r v as the virial radius of the halo, spherical collapse 
tells us that M = 47rr 3 A( 2 0/Qfc /3, wh ere A(z) is the over- 



density of collapse (see e.g. Henry 2000) and pb is the 
background matter density today. We use comoving co- 
ordinates throughout. By equating these two expressions, 
one can eliminate p s and describe the halo by its mass M 
and concentration c. 



Following Cooray et al (2000b), we take the concentra- 
tion of dark matter halos to be 



c(M,z) = a(z) 



M 



n -b(z) 



M*(z) 



(19) 



-0.3 



where a{z) = 10.3(1 + z)" - 3 and b(z) = 0.24(1 + z)~ 
Here M*(z) is the non- linear mass scale at which the peak- 
height threshold, u(M,z) = 1. The above concentration 
is chosen so that dark matter halos provide a reasonable 
match to the th e non-linear densi ty power spectrum as 
predicted by the Peacock fc Dodd£ (1996); it extends the 
treatment of EeijaE (2000) to the redshifts of interest for 
lensing. We caution the reader that eqn. ( p"9] ) is only a 
good fit for the ACDM model assumed. 

2.4. Results 

In Fig. [l](a-b), we show the density field power spectrum 
today (z = 0), written such that A 2 (fc) = k 3 P{k)/2w 2 
is the power per logarithmic interval in wavenumber. In 
Fig ^(a), we show individual contributions from the single 
and double halo terms and a comparison to the non-linear 



power spectrum as predicted by the Peacock & Dodd; 
(1996) fitting function. In Fig. Mb), we show the depen- 
dence of density field power as a function of maximum 
mass used in the calculation. 

Since the bispectrum generally scales as the square of 
the power spectrum, it is useful to define 



2 y/B{k,k,k), 



(20) 



which represents equilateral triangle configurations, and 
its ratio to the power spectrum 



Qcq{k) 



A c 2 q (fc) 
A 2 (fc) 



In second order perturbation theory, 



PT _ o _ 2 /63 

leq — x j iL m 



(21) 



(22) 



and under hyper-extend ed perturbation theory (HEPT; 
Scoccimarro fc Frieman 1999), 



C PT (fc) 



1 + 2"( fe ) +1 ' 



(23) 



which is claimed to be valid in the deeply nonlinear regime. 
Here, n(k) is the linear power spectral index at k. 



In Fig. yj(c-d), we show A 2 q (fc) separated into its various 
contributions (c) and as a function of maximum mass (d) . 
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Since the power spectra and equilateral bispectra share 
similar features, it is more instructive to examine Q cq (k) 
(see Fig. ^). Here we also compare it with the second or- 
der perturbation theory (PT) and the HEPT prediction. 
In the halo prescription, Q cq at k ;> 10fc nO niin ~ 10/iMpc -1 
arises mainly from the single halo term. The HEPT pre- 
diction exceeds the halo prediction on larger scales and 
falls short on smaller scales. 

2.5. Discussion 

Even though the dark matter halo formalism provides 
a physically motivated means for calculating the statistics 
of the dark matter density field, there are several limita- 
tions of the approach that should be borne in mind when 
interpreting the results. 

The approach assumes all halos to be spherical with a 
single profile shape. Any variations in the profile through 
halo mergers and resulting substructure can affect the 
power spectrum and higher order correlations. Also, real 
halos are not perfectly spherical which affects the config- 
uration dependence of the bispectrum. 

Furthermore, there are parameter degeneracies in the 
formalism that prevent a straightforward interpretation of 
observations in terms of halo properties. For example, 
one might think that the power spectrum and bispectrum 
can be used to measure any mean deviation from th e as 



sum ed NFW profile form. However as pointed out by Sel 



jak| (2000), changes in the slope of the inner profile can be 
compensated by changing the concentration as a function 
of mass; this degeneracy is also preserved in the bispec- 
trum. 

We do not expect these issues to affect our qualita- 
tive results. However, if this technique is to be used for 
precision studies of cosmological parameters, more work 
will be required in t esting it qu antitatively against sim- 
ulations. Studies by Ma fc Fry (2000a) and Scoccimarro 



et ajj (2000) show that the bispectrum predictions of the 



halo formalism are in good agreement with simulations, 
at least when averaged over configurations. The replace- 
ment of individual halos found in numerical simulations 



with synthetic smooth halos with NFW profiles by Ma & 



Fry| (2000b) show that the smooth profiles can regenerate 



the measured power spectrum and bispectrum in simula- 
tions. This agreement, at least at scales less than 10A: non ii n , 
suggests that mergers and substructures may not be im- 
portant at such scales. 

3. CONVERGENCE POWER SPECTRUM AND BISPECTRUM 

3.1. Power Spectrum and Variance 

The angular power spectrum of the convergence is de- 
fined in terms of the multipole moments K; m as 



( K *m K l'm'} = CfSu'Sr, 



(24) 



Ci is numerically equal to the flat-sky power spectrum in 
the flat sky li mit. It is related to the dark matter power 
spectrum by (|Kaiser| 1992; 1998) 



Cf = 



(25) 



where r is the comoving distance and dA is the angular 
diameter distance. When all background sources are at a 




10 u 10' 
k (h Mpc~ 

Fig. 2. — Qeq(k) at present broken into individual contribu- 
tions under the halo description and compared with second order 
perturbation theory (PT) and hyper-extended perturbation theory 
(HEPT). 



distance of r s , the weight function becomes 
W(r 



3 n Hi d A (r)d A (r s - r) 
2 m c 2 a d A (r s ) 



(26) 



for simplicity, we will assume r s = r(z s = 1) throughout. 
In derivi ng Eq. (p5[), wc have used the Limber approxi- 
mation ( Limber] 1954) by setting k = l/d A and the flat- 
sky approximation. In |Cooray et al (2000b), we used the 
projected mass of individual halos to construct the weak 
lensing power spectrum directly. The two approaches are 
essentially the same since the order in which the projection 
is taken does not matter. 

In Fig. ^, we show the convergence power spectrum of 
the dark matter halo s compared with that predicted by 
the iPeacock & Dodds (1996) power spectrum. The lensing 



power spectrum due to halos has the same behavior as the 
dark matter power spectrum. At large angles (I <; 100), 
the correlations between halos dominate. The transition 
from linear to non-linear is at I ~ 500 where halos of mass 
similar to Af*(z) contribute. The single halo contributions 

start dominating at / > 1000. 

As shown in Fig. ||(b), and discussed in Cooray et a\ 
(2000b), if there is a lack of massive halos in the observed 
fields convergence measurements will be biased low com- 
pared with the cosmic mean. The lack of massive halos 
affect the single halo contribution more than the halo-halo 
correlation term, thereby changing the shape of the to- 
tal power spectrum in addition to decreasing the overall 
amplitude. 

Similar statements apply to variance statistics (second 
moments) in real space. The variance of a map smoothed 
with a window is related to the power spectrum by 



(27) 



where Wi are the multipole moments (or Fourier transform 
in a flat-sky approximation) of the window. For simplicity, 
we will choose a window which is a two-dimensional top 
hat in real space with a window function in multipole space 
of Wi{a) = 2Ji(x)/x with x — la. 
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Fig. 3. — Weak lensing convergence spectra under the halo description, (a) Angular power spectrum broken into individual contributions, 
(b) Mass cut off effects on the angular power spectrum, (c) Equilateral bispectrum broken into individual contributions, (d) Mass cut off 
effects on the equilateral bisp<ectrum. Also shown in (a) is the prediction from the PD nonlinear power spectrum fitting function. The mass 
cuts are the same as in Fig. hi and we have assumed that all sources are at z a = 1. 



In Fig. |f|(a-b), we show the second moment as a function 
of smoothing scale a. Here, we have considered angular 
scales ranging from 5' to 90', which are likely to be probed 
by ongoing and upcoming weak lensing experiments. As 
shown, most of the contribution to the second moment 
comes from the double halo correlation term and is mildly 
affected by a mass cut off. 

3.2. Bispectrum and Skewness 
The angular bispectrum of the convergence is defined as 



{f^lxmi ^l^rai ^l^ra^ ) 



h h h 
mi TYI2 777,3 



(28) 



Exte nding our derivation of the SZ bispectrum in Cooray 
et al (2000b), we can write the angular bispectrum of the 
convergence as 



(2Zi + l)(2I 2 + l)(2I 3 + l) 



4tt 



h_ 

dA ' 



dA : 



h h h 


dA ' 



(29) 



The more familiar flat -sky bispectrum is simply the ex- 
pression in brackets ( |Hu| 2000). The basic properties 
of Wigncr-37 symbol introduced above can be found in 
Cooray et af (2000a). 

Similar to the density field bispectrum, we define 



(30) 



involving equilateral triangles in /-space. 

In Fig. |(c-d), we show A^ qi . The general behavior of 
the lensing bispectrum can be understood through the in- 
dividual contributions to the density field bispectrum: at 
small multipoles, the triple halo correlation term domi- 
nates, while at high multipoles, the single halo term dom- 
inates. The double halo term contributes at intermediate 
Vs corresponding to angular scales of a few tens of arcmin- 
utes. 

In Fig. |^, we show the configuration dependence 



K hh - 



2tt 



A 2 , 

cqi 



(31) 



as a function of h an d h when l 3 = 1000. The surface, 
and associated contour plot, shows the contribution to bis- 
pectrum from triangular configurations in I space relative 
to that from the equilateral configuration. Due to the tri- 
angular conditions associated with Vs, only the upper tri- 
angular region of h~h space contribute to the bispectrum. 
The symmetry about h — h line is due to the intrinsic 
symmetry associated with the bispectrum. Though the 
weak lensing bispectrum peaks for equilateral configura- 
tions, the configuration dependence is weak. 

The skewness is simply one, easily measured, aspect of 
the bispectrum. It is associated with the third moment of 
the smoothed map (c.f. eqn. p7y) 



1 

47T 



E 



(2i 1 + l)(2Z 2 + l)(2i 3 + l) 

47T 
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We then construct the skewness as 



(32) 



(33) 



The effect of the mass cut off is dramatic in the third 
moment. As shown in Fig |](c-d) , most of the contributions 
to the third moment come from the single halo term, with 
those involving halo correlations contributing significantly 
only at angular scales greater than ~ 25'. With a mass 
cut off, the total third moment decreases rapidly and is 
suppressed by more than three orders of magnitude when 
the maximum mass drops to 10 13 Mg. The skewness only 
saturates when the maximum mass is raised to a few times 
10 15 Mg. Even though a small change in the maximum 
mass does not greatly chang e the convergence power spec- 
trum (Fig. 3 of Cooray et al| 2000b), the third moment, or 
the bispectrum, is strongly sensitive to the rarest or most 
massive dark matter halos. 

In Fig. ^ we plot the skewness as a function of maximum 
mass, ranging from 10 11 to 10 16 Mg. Our total maximum 
skewness agrees with w hat is predicte d by numerical par- 
ticle mesh simulations ( White fc Hu| 1999) and yields a 
value of ~ 116 at 10'. However, it is low er than pr edicted 
by HEPT arguments and simulations of Jain et al (2000), 
which suggest a skewness of ~ 140 at angular scales of 
10'. The skewness based on second-order PT is factor of 
~ 2 lower than the maximum skewness predicted by halo 
calculation. As shown, the PT skewness decreases slightly 



from angular scales of few arcmins to 90' and increases 
thereafter. 

The effect of maximum mass on the skewness is inter- 
esting. When the maximum mass is decreased to 10 15 Mg 
from the maximum mass value where skewness saturates 
(~ 3x 10 15 Mg), the skewness decreases from ~ 116 to 98 
at an angular scale of 10', though the convergence power 




Fig. 5. — The bispectrum configuration dependence as a 

function of l\ and I2 with Z3 = 1000. Due to triangular conditions 
associated with Z's, only the upper triangular region in I1-I2 space 
contribute to the bispectrum. 
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Fig. 6. — The skewness, S:j(a), as a function of angular scale. 
ShoAjoi here is the skewness values with varying maximum mass as in 
Fig. [N. For comparison, we also sh ow sltpwnpss v alues as measured in 
particle-mesh (PM) simulations of White & Hu (1999), as predicted 
by hyper-extended perturbation theory (HEFT; dashed line) and 
second-order perturbation theory (PT; long-dashed line). 



spectrum only changes by less than few percent when the 
same change on the maximum mass used is made. When 
the maximum mass used in the calculation is 10 13 Mq, 
the skewness at 10' is ~ 8, which is roughly a factor of 15 
decrease in the skewness from the total. 

The variation in skewness as a function of angular scale 
is due to the individual contribution to the second and 
third moments. The increase in the skewness at angular 
scales less than ~ 30' is due to the single halo contribu- 
tions for the third moment. The triple halo correlation 
terms dominate angular scales greater than 50', leading to 
a slight increase toward large angles, e.g. from ~ 74 at 
40' to ~ 85 at 90'. However, this increase is not present 
when the maximum mass used in the calculation is less 
than ~ 10 14 Mq. Even though mass cut off affects the 
single halo contributions more than the halo contribution, 
at such masses, the change in halo contribution with mass 
cut off prevents an increase in skewness at large angular 
scales. 

The absence of rare and massive halos in observed fields 
will certainly bias the skewness measurement from the cos- 
mological mean. One therefore needs to exercise caution 
in u sing the s kewness to co nstrain cosmological models 
" 1999). 



Hui 



In Cooray et al 



(2000b), we suggested that 
lcnsing observations in a field of ~ 30 deg 2 may be ad- 
equate for an unbiased measurement of the convergence 
power spectrum. For the skewness, observations within a 
similar area may be biased by as much as ~ 25%. This 
is consistent with the sampling errors found in numerical 
simulations: l a errors of 24% at 10' with a 36 deg 2 field 
(White & Hu 1999). To obtain the skewness within few 
percent of the total, one requires a fair sample of halos 
out to ~ 3 x 10 15 Mq, requiring observations of ~ 1000 
deg 2 , which is within the reach of upcoming lensing sur- 
veys involving wide-field cameras, suc h as the MEGA CAM 
at Canada- France-Hawaii- Telescope ( Bouladc ct aj 1998), 
and proposed dedicated telescopes (e.g., Dark Matter Tele- 
scope; Tyson, private communication). 



Still, this does not mean that non-Gaussianity measured 
in smaller fields will be useless. With this halo approach 
one can calculate the expected skewness if one knows that 
the most massive halos are not present in the observed 
fields. This knowledge may come from external informa- 
tion such as X-ray data and Sunyaev-Zel'dovich measure- 
ments or internally from the lensing data. 

3.3. Related Statistics 

The halo description in general allows one to test the 
effect of rare massive halos on any statistic related to the 
two and three point functions. In particular, it can be 
used to design more robust statistics. 

Generalized three poin t statistics have been considered 
previously by Jain et al| (200 0) following Nusser fc Dckc] 
(1993) and Juszkiewicz ct a] (1995). One such statistic is 



the (k|k|) k>0 , which is expected to reduc e the sam pling 
variance from rare and massive halos (see, Jain et a] 2000 

1/2 

for details). This statistic is proportional to ^« 3 ) / ^k 2 ) . 
In Fig. |^(a), we show this statistic as a function of max- 
imum mass used in the calculation. We still find strong 
variations with changes to the maximum mass. Similar 
va riations we re also present in other statistics considered 
by |Jain et al| (2000). 




35 45 55 65 

a (arcmin) 

Fig. 7.t— Generalized skewness statistic S-j(cr, m). (a) m = 1/2 
following Jain et al (2000). (b) m = 3.7 chosen to minimize the 



mass cut off dependence. 
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Consider instead the generalized statistic 



(34) 



where to is an arbitary index. We varied m such that 
the effect of mass cuts are minimized on skewness. In 
Fig. |?](b), we show such an example with to = 3.7. Here, 
the values are separated to two groups involving with most 
massive and rarest halos and another with halos of masses 
10 14 Mq or less. Though the values from the two groups 
agree with each other on small angular scales, they depart 
significantly above 25' reaching a difference of 2.5 at 80'. 
Statistics involving such a high index to, weight the single 
halo contributions highly when the most massive halos are 
present, whereas they weight the halo correlation terms 
more strongly for M < 10 14 Mq. To some extent this 
may be useful to identify the presence of rare halos in the 
observations. 

However the consequence of using these generalized 
statistics is that one progressively loses their independence 
on the details of the cosmological model, e.g. the shape 
and amplitude of the underlying density power spectrum, 
as one departs from m = 2, thereby contaminating the 
probe of dark matter and dark energy. Further work is 
necessary find the optimal trade off between robustness 
and cosmological independence of these and other gener- 
alized statistics. 

Another observable statistic is the probability distribu- 
tion function (pdf) of the convergence maps smoothed on 
the scale a. This poss ibility has been recently studied by 
Jain & van Waerbeke (1999), where the reconstruction of 
pdf using peak statistics were considered. Using the Edge- 
worth expansion to capture small deviations from Gaus- 
sianity, one can write the pdf of convergence to second 
order as 



P(k) 



1 



1 



3 - K ( CT ) 2 /2(k 2 (<t)) 



(35) 
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Fig. 8. — The probability distribution function of the weak lensing 
convergence as a function of maximum mass used in the calculation 
at an angular scale of 12'. From top to bottom, the curves range 
from 10 11 to 10 16 M . 



where H3 (x) 
mial (see, 
In Fig. 



3x is th e third order Hermite polyno- 



Juszkiewicz et aj 1995 for details). 
SI, we show the pdf of convergence at 12' as a 
function of maximum mass used in the calculation. As 
shown, the greatest departures from Gaussianity begin to 
occur when the maximum mass included is greater than 
10 14 Mq. Given that we have only constructed the pdf 
using terms out to skewness, the presented pdfs should 
only be considered as approximate. With increasing non- 
Gaussian behavior, the approximated pdfs are likely to 
depart from this form especially in the tails. As studied in 
Jain & van Waerbeke (1999), the measurement of the full 
pdf can potentially be used a a probe of cosmology. Its 
low order properties describe deviations from Gaussianity 
near the peak as opposed to the skewness which is more 
weighted to the tails. 

4. SUMMARY & CONCLUSIONS 

We have presented an efficient method to calculate the 
non-Gaussian statistics of lensing convergence at the three 
point level based on a description of the underlying den- 
sity field in terms of dark matter halos. The bispectrum 
contains all of the three point information, including the 
skewness. The prior attempts at calculating lensing bis- 
pectrum and skewness were limited by the accuracy of 
perturbative approximations and the dynamic range and 
sample variance of simulations. 

Though the present technique provides a clear and an 
efficient method to calculate the statistics of the conver- 
gence field, it has its own shortcomings. Halos are not 
all spherical, which can to some extent affect the configu- 
ration dependence in moments higher than the two point 
level. Substructures due to mergers of halos can also intro- 
duce scatter. Though such effects unlikely to dominate our 
calculations, further work using numerical simulations will 
be necessary to determine to what extent present method 
can be used as a precise tool to study the higher order 
statistics associated with weak gravitational lensing. 

The dark matter halo approach also allows one to study 
possible selection effects that may be present in weak lens- 
ing observations due to the presence or absence of rare 
massive halos in the small fields that are observed. We 
have shown that the weak lensing skewness is mostly due 
to the most massive and rarest dark matter halos in the 
universe. The effect of such halos is stronger at the three 
point level than the two point level. The absence of mas- 
sive halos, with masses greater than 10 14 Mq, leads to 
a strong decrease in skewness, suggesting that a straight- 
forward use of measured skewness values as a test of cos- 
mological models may not be appropriate unless prior ob- 
servations are available on the distribution of masses in 
observed lensing fields. 

One can correct for such biases using the halo approach, 
however. To implement such a correction in practice, fur- 
ther work will be needed to calibrate the technique pre- 
cisely against simulations across a wide range of cosmolo- 
gies. Efficient techniques to correct for mass biases both in 
the lensing power spectrum and bispectrum will be needed. 
Alternatively, this technique can be used to search for gen- 
eralized three point statistics that are more robust to sam- 
pling issues. Given the great potential to study the dark 
matter distribution through weak lensing, this issues merit 
further study. 
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